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ABSTRACT 

We present the results of high resolution spectroscopy of the prototype polar 
AM Herculis observed with Chandra High Energy Transmission Grating. The 
' X-ray spectrum contains hydrogen-like and helium-like lines of Fe, S, Si, Mg, Ne 

O 



> 
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and O with several Fe L-shell emission lines. The forbidden lines in the spectrum 
are generally weak whereas the hydrogen-like lines are stronger suggesting that 
emission from a multi-temperature, collisionally ionized plasma dominates. The 
helium-like line flux ratios yield a plasma temperature of 2 MK and a plasma den- 
sity 1-9 xlO 12 cm -3 , whereas the line flux ratio of Fe XXVI to Fe XXV gives an 
ionization temperature of 12. keV. We present the differential emission mea- 
sure distribution of AM Her whose shape is consistent with the volume emission 



measure obtained by multi-temperature APEC model. The multi-temperature 
plasma model fit to the average X-ray spectrum indicates the mass of the white 
dwarf to be ~ 1.15 M . 

From phase resolved spectroscopy, we find the line centers of Mg XII, S XVI, 
resonance line of Fe XXV, and Fe XXVI emission modulated by a few hundred 
to 1000 km s _1 from the theoretically expected values indicating bulk motion of 
ionized matter in the accretion column of AM Her. The observed velocities of 
Fe XXVI ions are close to the expected shock velocity for a 0.6 M white dwarf. 
The observed velocity modulation is consistent with that expected from a single 
pole accreting binary system. 



Subject headings: novae, cataclysmic variables - stars:individual-(AM Her) - X- 
rays: binaries 
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Introduction 



Polars are a class of magnetic cataclysmic variables that are synchronously locked, close 
interacting binary systems with a highly magnetized white dwarf primary. The high magnetic 
field that can range from 10 MG to 230 MG prevents the formation of an accretion disk and 
channels the accreting matter along the field lines to the surface of the white dwarf. The 
slowing down of supersonic accreting matter forms a shock near the surface of the white 
dwarf. The shock heats up the in-falling matter to very high temperatures (~ 10 8 K) 
and the hot plasma in the post shock region cools off by emitting hard X-rays via thermal 
bremsstrahlung process. A part of these hard X-rays heat up the surface of the white dwarf 
which in turn emits soft X-rays. 



AM H er, is the pro t otype of polars at a distance of 91 pc (IGansicke et al.lll995l ). Ac- 
cording to lYoung et al.l (ll98ll ) AM Her has an orbital period of 3.1 hours and consists of 
a white dwarf of mass 0.39 M Q and radius 0.26 R & , and an M dwarf secondary with a 
mass of O.26M . However, various other estimates for the mass of the white dwarf reported 
in the literature quote values in the range of 0.5 M^ to 1.22 Mm (IGansicke et al. 



Wu. Changmugam fc Shavivll 19951 ; iMukai fc Charleslll987l ; iMouchet 
1998h. 
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1998 



Cropper. Ramsay &: Wu 



High resolution X-ray spectroscopy with Chandra 's High Energy Transmission Grating 
(HETG) allows us to resolve many spectral lines which in turn allow us to study the X- 
ray emission processes, and temperature, density and ionization properties of the X-ray 
plasma. The high resolving power of Chandra can also help us to determine the Doppler 
shifts of emission lines, which can be used to study the kinematics of X-ray emitting region 
in Cataclysmic Variables. 

Here, we present the results of high resolution X-ray sp ectroscopy of AM Her based 
on the analysis of archival Chandra data taken with HETG (ICanizares et al.ll2005l ). The 
organization of the paper is as follows. In SJ2], we describe the observations, data analysis 
and the search for the presence of Doppler shifts of different emission lines, discussion in §21 
and conclusions in £ffl 



2. Observations and Analysis 

Chandra observed AM Her on 2003 August 15 (Obsld 3769) using HETG for 94 ks. The 
data were continuous and covered ~8.5 orbital cycles. The HETG obs ervations were taken 
in combination with the Advanced CCD Imaging Spectrometer (ACIS; iGarmire et al.l 120031 ) 
in the faint mode. The HETG carries two transmission gratings: Medium Energy Grating 
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(MEG) and High Energy Grating (HEG). The absolute and relative wavelength accuracy of 
HEG are ±0.006 A and ± 0.001 A respectively (Chandra proposers' observatory guide v. 7). 
This observation was previously used in a comparative study of iron Ka complex in five 
of th e magnetic cataclysmic variables observed by Chandra at that time (IHellier fe Mukai 
2004T ). 



AM Her generally remains in a high state (V m ~ 13 mag), dropping into a low state 
(V m > 15 mag) occasionally. Th e observations repo rted by the American Association of 
Variable star observers (AAVSO; iMattei et al.l Il980l ) near the Chandra observation date 
show that the star was in an intermediate state during the HETG observations (V m ~ 14). 

We used the Chandra Interactive Analysis of Observations (CIAO) v3. 2 software pack age 
for data reduction and the extraction of the spectra, and XSPEC vl2.2 (jArnaudlll996l ) for 
spectral analysis. The lightcurves extracted from un-dispersed zeroth order HETG image 
show eclipses. There is no X-ray flaring activity in the source or in the background during 
the times of observations, and hence, no data were rejected. 

Using the tg extract tool available in CIAO, the spectra for -1 and +1 orders of HEG 
and MEG were extracted and summed together using add_grating_orders script to improve 
the statistics. We used the default binning of 0.0025 A for HEG and 0.005 A for MEG. 
Background subtraction is not applied to the data, since the background counts are less 
than 1% of the source counts. The average Chandra HEG and MEG spectra of AM Her 
are plotted in Figure [H and Figure [2] respectively. For clarit y, the data were smo othed by 
a Gaussian of width 0.0025 A (HEG) and 0.005 A (MEG) fauenemoerderl l2003l l Several 
emission lines are seen, an d were identified using Astrophysical Plasma Emission Database 
(APED; ISmith et al1l200lh . The lines identified are hydrogen-like lines of Fe XXVI, S XVI, 
Si XIV, Mg XII, Ne X and O VIII, helium-like triplets of Fe XXV, Ne IX and O VII, fluorescent 
Fe I line and several Fe L-shell lines. The hydrogen-like lines of S, Si, Mg, Ne and O are 
prominent in the spectrum while their corresponding helium-like lines are comparatively 
weaker or nearly absent altogether. The forbidden lines in hydrogen-like triplets of different 
ions are very weak compared to the corresponding resonance and inter-combination lines. 
Also, the hydrogen-like and helium-like lines of Fe are very strong along with a Fe fluorescent 
line signifying the presence of both high temperature plasma and photoionization in the 
spectrum. 
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2.1. Global spectrum analysis 



The AM Her spectrum is well fitted with four component APEC models, with a neutral 
absorber and a partial covering absorber with a covering fraction of 62%. The best fit values 
of temperatures from the multi-temperature APEC global fitting and their corresponding 
volume emission measures are summarized in Table [TJ The multi-temperature APEC models 
give a reasonable fit to the spectrum, providing a discrete representation of the temperature 
distribution of the plasma. In reality a continuous distribution of temperature is expected 
representing the physical conditions that exist in the X-ray emitting regions of magnetic 
cataclysmic variables (MCVs). Therefore, we have fitted multi-temperature plasma model 
to AM Her spectrum, as described below. 



The multi-temperature plasma model uses the calculations of lAizul (119731 ) to predict 
the temperature and the density of a magnetically confined hot plasma in the post-shock 
region. Full formulation and assum ptions used in the multi-temperature plasma model is 
described by ICropper et al.l (119991 ) . The multi-temperature plasma model was obtained 
from G. Ramsay via private communication and added to the XSPEC as a local model. The 
multi-temperature plasma model and a partial covering absorption component along with 
an absorption due to interstellar matter are used to fit the average spectrum of AM Her. 
To model the fluorescent Fe line a Gaussian centered at 6.4 keV is included. The ratio 
of cyclotron to bremsstrahlung cooling (eo) is fixed at 20 corresponding to the observed 
magnetic field (~14 - 22 MG) of AM Her. Varying the eo value over 10 - 30 does not affect 
the other spectral par ameters significantly. The mass accretion rate is fixed to 0.4 g s _1 cm -2 
( IGansicke et al.lll995l ) and the radius of the accretion column is fixed to 10 8 cm, typical of 
white dwarves. Elemental abundances are assumed to be solar, and the accretion column is 
divided into 100 vertical grids with a viewing angle of 35°. The best fit value for the mass 
of the white dwarf is found to be 1.15 ± 0.05 M Q with a reduced Xmin °f 0-5 for 2875 dof. 
The mass obta ined by using multi-temperatur e plasma model is consistent with the mass 
determined by ICropper. Ramsay fe Wul (119981 ). However, in this model authors assumed 
a one-temperature flow in post shock region of MCVs in w hich the electrons and the ions 
share the same temperatures. Recently, ISaxton et al.l (120051 ) improved on this and performed 
two-temperature hydrodynamic calculations for the post shock flow in accretion column of 
MCVs. Two-temperature flow is exp ected to occur in MC Vs with strongly magnetized (>20 



MG) massive white dwarf (~1 M Q ). ISaxton et al.l (120051 ) have compared the results of two- 



temperature calculations with the one-temperature calculations and found that the white 
dwarf mass estimated by using two-temperature hydrodynamic spectral models are lower 
than those obtained by one-temperature models. This bias is found to be more severe for 
massive white dwarf systems. 
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2.2. Line-by-Line Analysis 

In another method, we used Gaussians to fit the individual emission lines with the 
aim of measuring the emission line properties. This method has the virtue of not assuming 
plasma parameters inherent in MEKAL or APEC models, but to derive the plasma properties 
from the ratios of the fluxes estimated by using Gaussians. Both HEG and MEG spectra 
were fitted simultaneously. The continuum was estimated by an unabsorbed bremsstrahlung 
component of temperature 15.5 keV. Both the continuum and Gaussians were folded with 
the instrumental response. The normalization of the bremsstrahlung component was frozen 
after fitting the continuum using line free regions. The width of each line was tested for 
significant non-zero value by initially varying it for each line. It was found that the line 
widths, in general, were consistent with the resolution of the instrument. Therefore, the 
widths were fixed at zero to mimic the unresolved lines, except for the fluorescent line of 
Fe where the width was allowed to vary. The line center energy and the photon flux were 
allowed to vary for each line during the final fitting except for few weak lines for which the 
line-centers were frozen at the respective theoretical energies. Cash statistics was used to 
estimate the parameter values and confidence ranges. Though Cash statistics is a better 
criteria to determine the confidence range on a best fit parameter when the data bins have 
few counts, in its original form, Cash statistics could not provide a goodness of fit similar to 
X 2 ■ XSPEC uses a modified function which provides a goodness-of-fit criterion similar to x 2 
in the limit of larger counts. The background counts are not subtracted from source counts 
while modelling, as the background counts form only 1% of the source spectrum counts. Not 
subtracting the background counts from the source counts does not affect our interpretation. 
The best fit model gives a C-statistic value of 7747 for 7875 dof. Table [2] lists the fluxes and 
the equivalent widths of emission lines along with their probable identifications. The model 
spectrum is plotted as a dashed line in Figure [1] and Figure [2j The individual emission lines 
are also marked in the spectrum. 

In the average spectrum, the line-centroids of several lines show a significant deviation 
from the predicted wavelength of the lines. Assuming the observed shifts are real and not 
an artefact, we interpret the shifts as the shifts caused due to Doppler effect. A comparison 
of the observed line centers with the theoretically predicted line centers thus yields different 
Doppler velocities for Fe XXVI, resonance line of Fe XXV, These Doppler shifts are listed in 
the second column of Table 
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2.3. Emission measure analysis of the average spectrum 



The volume emission measure (VEM) is the measure of the "amount of material" avail- 
able in a plasma to produce the observed flux, which gives us an idea of how the emitting 
material is distributed with the emitting temperature. For an optically thin plasma in col- 
lisional ionization equilibrium, the relation for VEM can be written as ([Griffiths fe Jordan 
19981 ). 



n e n h dV 



[max) 



,2 fi 



where, fi is observed flux in a line feature, Gi is emissivity {photons cm 3 s 1 ), d is distance 
to the source, n e and are the electron and hydrogen densities in cm -3 , respectively. 

Using Equation [TJ we can calculate the maximum emission measure as a function of 
temperature for each emission line. The combination of emission measure of individual lines 
can be used to constrain the emission distribution of the source. The VEM estimated using 
Equation [J f or different lines a re plotted in Figured The emissivities of the lines are taken 
from APED (ISmith et al.ll200ll ). The curves are the loci corresponding to the VEM required 
to produce the observed flux from an isothermal plasma as a function of plasma temperature. 
At a given temperature, a point on one of these curves represents a maximum EM that the 
plasma can have at that temperature. The solid lines in the figure represent hydrogen-like 
lines of 0, Ne, Mg, Si, S and Fe, dashed lines correspond to resonance lines of VII, Ne IX 
and Fe XXV, whereas the dash-dot-dash lines represent the VEM of Fe L-shell lines. 

The VEM defined in EqnJTJ in the logarithmic differential form is called as Differential 
Emission Measure (DEM). The DEM gives a correlation between the amount of emitting 
power and the amount of emitting material in the plasma as a function of temperature. The 
DEM was estimated by performing a M arkov Chain Monte-Carlo anal ysis using a Me tropolis 
algorithm (MCMC[M]) JKashvaplll998h implemented in PINTofALE faashvapll^OoT ) on the 
set of sixteen brightest lines in the AM Her spectrum. The MCMC(M) method gives an 
estimate of emission measure distribution over a pre-selected temperature region with the 
DEM defined for each bin. Here, we used a temperature grid ranging from logT = 6.0 
to logT = 8.8 with AlogT = 0.2. The lower and the upper limit of the temperature 
region are choosen to represent the hydrogen and helium-like lines of O and, hydrogen and 
helium-like lines of Fe respectively, present in the spectrum. The intermediate temperatures 
are constrained by the L-shell lines of Fe, and helium-like lines of Ne, Mg, S and Si. The 
reconstructed DEM is plotted in Figure[3]with 95% confidence limits shown as shaded region. 
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2.4. Phase dependent variation of spectral lines 



To look for modulations in the continuum and emission lines with the binary period, the 
data were f olded u s ing th e eph emeris for the magn etic phase and a period of 0.128927 days 
as given by iTapial (Il977l ). and lYoung et al.l (Il98ll ). Two methods were adopted to obtain 
the phase resolved spectra. In the first method, the spectra were extracted for five non- 
overlapping regions of 0.2 phase interval. Except for the interval centered around the orbital 
phase minima, other intervals had enough number of counts to estimate the line energy and 
the photon flux in many of the bright emission lines. From these spectra it was observed 
that while some lines showed shifts in the line energy with phase, others did not. The 
second method is based on extracting the spectra into continuous overlapping regions (See 
Hoggenwerf et al.l 12004 ) . In this method, the spectra were accumulated for twenty phase 
intervals of width A$ = 0.25, with phase intervals starting at $o — 0.0, 0.05, 0.1 . . ., 
where, $ is the orbital phase. Spectra obtained from both the methods were analysed 
for their continuum and emission line properties. The continuum was modelled using a 
bremsstrahlung after masking the regions of emission lines identified using APED. 

Different regions of spectra were examined for studying the line emission properties 
as a function of phase. Figure H] shows spectra for four non-overlapping phase bins in the 
wavelength region 1.75 - 2.0 A which covers the Fe XXVI line, Fe XXV (r, i, and /), and 
the fluorescent line of Fe I. The best fit line energies were compared with the reference 
line energies to determine the shifts of different emission lines observed. The Fe XXVI line 
consists of a doublet, Lycti and Ly«2 at 6.973 and 6.952 keV respectively. The two lines can 
not be resolved with Chandra HEG resolution of ~ 35eV at 7 keV. The branching ratio 2:1 
of th e Ly«! and Lya2 yield a line centroid of 6.966 keV (IPike et al.lll996l ; iHellier fe Mukai 
2004 ) . This value was used as a reference to determine any shift in the Fe XXVI line energy. 
The stabilizing transitions of doubly excited ions lead to dielectronic recombination satellite 
(DES) lines redward of Fe XXVI lines. The Fe XXVI lines are mainly affected by a feature 
at 6.92 keV whos e intensity is ~ 8% of the Fe XXVI lines at 5 x 10 7 K and falls below 
5% at 8 x 10 7 K JPubau et al.lll98lh . Hence the contribution of DES lines to Fe XXVI in 
AM Her can be neglected due to high plasma temperature of 1.8 x 10 8 K (~15.6l 3 *j keV; 
the maximum temperature from 4T APEC fit from Tabled]). The intensities of DES lines 
depend both on temperature and intensity of adjacent ionization states, specially Fe XXIV 
and Fe XXIII. At temperatu res above 3.5 x 10 7 K, the contribution of DES is negligible 
(jOelgoetzty fe Pradhanll200ll ). and hence, we did not consider the effects of DES lines. To 
include the contribution of unresolved % lines and / line, we used three Gaussians representing 
the r, % and / lines. For the calculation of line shifts in Fe XXV, we used the rest frame 
energy of r at 6.700 keV, % at 6.675 keV and / at 6.636 keV. 
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The strength of DES lines approximately scales as Z A and hence, their effects have been 
neglected in the analysis of S, Si and Mg. We used 2.623, 2.005 and 1.473 keV as reference 
line energies for S XVI, Si XIV and Mg XII respectively. 

From the phase resolved spectra, we found that the line centers of individual lines for 
different phase bins deviate from the predicted wavelengths. The difference between the 
observed and expected emission line center is used to calculate the Doppler shift of the 
emission line from different phase bins. 

The calculated Doppler shifts of emission lines of Fe XXVI and Fe XXV (r) as a function 
orbital phase are plotted in Figure [5j Doppler shifts obtained from non-overlapping phase 
bins are shown as filled circles. The shifts from overlapping phase bins are shown for a better 
visualization only. The figure clearly shows a non-zero shift in the two lines. The amplitude 
of Doppler shift of different lines also shows a variation with orbital phase. The modulation 
in the shift in line center with phase is most clearly seen in Fe XXVI line as compared to 
the Fe XXV line. In order to quantify the observed modulation of the Doppler shifts with 
phase in the above mentioned lines a sinusoid+constant model was fitted to the data using 
X 2 minimization. To calculate the un-modulated velocity shift, we fit a sinusoid with non- 
zero mean to the observed modulation of all lines using non-linear least square fitting. We 
have also checked for the presence of any Doppler shifts in fluorescent Fe, S XVI and Mg XII 
lines. The mean Doppler velocity determined from the average spectrum of different lines are 
summarized in the second column of Table [3J The third and the fourth columns of Table [3] 
list a constant velocity shift and the amplitude of the sinusoid, respectively. The sinusoid 
best fit for different lines are superposed over the doppler shifts as dashed lines in Figure [5j 
For Fe XXVI we get a semi-amplitude of 790 ± 40 km s" 1 with a mean of 220 ±26 km s _1 . 
Just to compare, we fit a constant model to the velocity shifts and list them in column five 
of Table [31 A constant fit to the Fe XXVI data gives a velocity shift 565±195 km s _1 , but 
the fit is considerably worse with a xl — 3.29 as compared to a xl — 0.02 for the sinusoid fit. 
The shifts in the Fe XXV (r), S XVI and Mg XII are however well fit with a constant velocity 
shift of 770 ± 75 km s^ 1 (xl = 0.06), 500±160 km s" 1 (xl = 0.79), and 280±195 km s" 1 
(xl = 0.15) respectively. The xl values for the constant fits and sinusoid fits are similar. 
The fluorescent Fe line shows no detectable shift, with 2cr upper limit of 400 km s _1 . The 
Si XIV line is too weak to determine the parameters accurately, therefore, no constraints are 
put on the shift of this line. 

Although no modulation is apparent in the velocity shifts of the Fe XXV (r), S XVI, 
and Mg XII, the error bars on the shifts of their lines can not rule out the possibility of 
modulation being present at some level. To estimate the level of modulation, we have force 
fitted a constant plus a sinusoid component similar to the one used for Fe XXVI. The semi- 
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amplitude of the sinusoid and the constant value for the velocity shifts of Fe XXV, S XVI 
and Mg XII are also listed in Table [3j The amplitude of modulation in these lines are of the 
order of the error bars. 



Discussion 



3.1. Average X-ray spectrum 

Ratios of the line fluxes of helium-l ike lines provide good diagnostics of density an d 
temperature of the line forming regions (IGabriel fe Jordan! Il969l ; iPorquet fe Dubaull2000l ). 
The transition from the excited 1 Pi, 3 Pi and 3 Si levels to the ground level x Pi forms the three 
strongest lines of helium triplets: the resonance line (r), the intercombination line (i) and the 
forbidden line (f) respectively. The analytical relations between th e electron density R(n P ) 
and temperature G(T P ) axe, R(n e ) = f/i and G{T e ) = (f + i)/r (IGabriel fe Jordan! Il969 
Bluementhal et all 11972 ) . 



The He-like triplets of oxygen and neon are used as the diagnostics of low temperature 
regions. As we could not resolve the He-like lines of Fe, we do not use Fe XXV as a plasma 
diagnostic. We can only put an upper limit on the fluxes of 'f ' lines of Ne IX and O VII and 
hence, we could only estimate the upper limits of the two ratios. The helium-like lines of Si, 
S and Mg are too weak in the spectrum and hence do not allow us to put any constraints 
on their G and R ratios. Using the values listed in Table [2J we get R < 0.48 and G < 0.76 
for oxygen triplet and R < 0.5 and G < 0.8 for neon triplet, as 3<r upper limits. 

Comparing the measured G(T e ) values of O VII and Ne IX with the theoretical relation 
between G(T e ) and electronic temperatures (IPorquet fe Dubaull2000l ) suggests a temperature 
greater than ~ 2MK. This temperature agrees with the lowest temperature of 0.13 keV 
(1.5 MK) of the four APEC components. Similarly comparing the measured values of R(n e ) 
with theoretically predicted values for a collisional dominated hybrid plasma implies a density 
greater than 2 x 10 12 cm" 3 . Since the temperature and density of the O and Ne helium-like 
line emitting regions are very close, we assume that the emission of the two triplets takes 
place very close to each other. However, it should be noted that the presence of UV radiation 
fields can also mimic high densities as the transition from / to i ca n also be triggered by 



UV photons. Stron g UV radiation is known to be present in AM Her (IWesemael et al.lll980 
Szkodv et al.1ll982h . 



The measured ratio of H- to He-like line intensities can be used to constrain the ioniza- 
tion temperature of the emitting plasma. We have used Fe ions to constrain the maximum 
ionization temperature in the post-shock region. A value of I(Fe XXVI)/I(Fe XXV) = 1.22 
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+0 ;ll suggests an ionization temperature of 12. keV (IMewe et allll985l ). This matches 



-0.: 



well with the maximum electron temperature 15. 6tJ]i keV of AM Her obtained from 4-T 
APEC global fit. The agreement with the continuum temperature and ionization temper- 
ature is consistent with the results obtained using moderate energy resolution ASCA data 



fllshida et al.lll997h . 



Volume emission measure analysis using individual lines shows that the plasma has a 
range of temperature from 1.5-100 MK. The absence of the N VI He-like triplet indicates an 
VEM distribution dominated by high temperatures. We obtain the lowest VEM from Fe XVII 
with 2.96xl0 50 cm" 3 at logT = 6.8 (~6 MK), whereas the peak VEM of 9.73xl0 53 cm" 3 
from Fe XXVI is at logT = 8.2 (~63 MK). We also estimated the differential emission 
measure using the fluxes of sixteen strong lines. We obtain a well constrained DEM(T) 
between logT = 6.0 and logT = 8.2. The minimum of the DEM is around logT = 6.6 — 6.8 
and there is an indication of a DEM maximum around logT = 8.2, but, the lack of sufficient 
number of lines in this region makes it ambiguous. Though the trial DEM ranges from 
logT = 6.0 to logT = 8.8, the lack of spectral information beyond logT = 8.2 makes it 
un-realistic. The 4-T APEC model provides only an approximate discrete representation of 
the continuous VEM distribution. The temperature distribution thus obtained is consistent 
with that determined using individual lines. The temperature with the lowest VEM/DEM 
determined by the former method agrees with the lowest VEM of 4-T model. 



3.2. Phase resolved spectroscopy 

According to the standard models of accretion in polars, supersonic matter in an accre- 
tion column forms a shock near the surface of the white dwarf heating up the post shock 
region to very high temperatures (10 8 K) thus emitting hard X-rays. As the post shock mat- 
ter falls toward the white dwarf surface, it cools down while getting slowed down. A phase 
dependent change in the spectral shape of the continuum is consistent with a non-uniform 
temperature distribution expected in the accretion column of polars. The central energies of 
emission lines of different elements, particularly hydrogen-like ions of Mg, S and Fe appear 
to be shifted from their theoretically expected values, implying possible bulk motion in the 
post-shock region of AM Her. The velocity shift increases as the atomic number increases 
from Mg to Fe. The observed gradient in the velocities of different ionization states of dif- 
ferent elements can be interpreted as a signature of cooling and the slowing down of the 
matter. 

The physical conditions in the accretion column can be described by velocity (v), tem- 
perature (T), density (n e ), height and size of the accretion column (h, r). At the shock 
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region, the values of th ese parame ters depend only on the mass, radius and the accretion 
rate of the white dwarf (jAizulll973l ). The velocity shift observed in emission lines of different 
ionization states of different elements seen in the Chandra HEG spectra of AM Her gives 
us a unique tool for the determination of the structure of the accretio n column of AM Her, 
using the relationships among the above listed parameters as given by lAizul (119731 ). 



The accretion column parameters in the shock region, calculated using different mass 
estimates of AM Her ranging fr om 0.39-0.91 M rep orted in the literature, and an accretion 



rate of m = 0.4 g cm 2 s 1 ( iGansicke et al.l Il995l ) are listed in Table HI The radius of 



the whi te dwarf for each mass value was calculated using the mass-radius relation for white 
dwarfs (jNauenbergill972l ) . and found to be in the range of 0.016-0.009 R®. Using these values 
we estimated the values of parameters like shock veloci ty, v. q h, shock tempe rature, kT sh and 
shock height, h, as listed in Table H] (See Appendix A, iTerada et al.l 1200 ll . for the relations 
of shock parameters with white dwarf mass, radius and accretion rate). We assumed the 
solar abundances for the elements and mean molecular weight of the in-falling gas, /i = 
0.615, which corresponds to a mixture of hydrogen and helium in the ratio of 0.7 and 0.3 
respectively. The calculated v s h for < 0.6 M is much less than the observed velocity 
shift of the Fe XXVI line. Thus, assuming the observed velocity of Fe XXVI as equal to v s h, 
we obtain the mass of the white dwarf to be ~ 0.6 M . This also puts some constraints 
on the geometry of the accretion column of AM Her by providing corresponding values of 
height and radius of the shock region. For a 0.6 M white dwarf, we estimated the heights at 
which different emission lines are emitted in the accretion column and list them in the sixth 
column of Table [3J Figure [6] shows the schematic of the distribution of different line emitting 
regions in the accretion column from the white dwarf surface. Different shades/patterns 
in the figure correspond to different emission lines labelled in the figure and the width of 
each region corresponds to 90% confidence levels of the height of the emitting region for the 
particular ion. 

A rough estimate about the upper limit on the white dwarf mass can be obtained using 
the fact that the heights at which diff erent emission lines are em i tted i n the accretion column 
depend on the mass of white dwarf. IWu. Cropper &: Ramsayl (120011 ) have investigated the 
ionization structure of the post shock region of MCVs and predicted line emissivity profiles for 
Fe XXVI and Fe XXV lines for white dwarf masses of 1.0 M and 0.5 M assuming different 
magnetic field strengths (eo = 0,1,10,100). They show that highly ionized Fe lines are 
emitted close to the shock front for a low mass white dwarf with high magnetic fields, whereas 
they are formed well below the shock front for high mass white dwarfs. Though the emissivity 
structure depends on the white dwarf mass, the ratio of these heights will be independent of 
the mass of the white dwarf. Based on the observed velocities for Fe XXVI and Fe XXV line 
emitting regions, we calculate the ratio of the heights of these two lines and compare this 
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with that obtained from line emissivity profiles (Figures 6 and 7 of IWu. Cropper fe Ramsay 



200fl ). Here we assume that the height of the maximu m emissivity region of a line coin cides 



with the height determined from its bulk velocity. From lWu. Cropper fe Ramsayl (120011 ). the 
ratios are found to be ~ 3 and ~ 5 for white dwarf of mass 1.0 M and 0.5 M respectively. 
The ratio of the heights of these two lines in AM Her is ~ 2 indicating a mass closer to 1 M Q . 

Using the mean emissivity over the region of maximum emissivity of the four emission 
lines Fe XXVI, Fe XXV(r), S XVI and Mg XII taken from APED and adopting a distance of 
91 pc for AM Her, we obtain an average VEM of 2.54xl0 54 cm -3 . Using this VEM and the 
values of kT s h and h listed in Table HI we determine the electron density, n s e h at the shock 
and the radius of the shock region, r^, and list them in Table HI 

In polars, the accreting matter is assumed to be streaming along the magnetic poles. 
The bulk motion velocity (i>) of this matter will be seen from different angles along the line 
of sight due to the rotation of the white dwarf. The apparent velocity of accreting matter 
at any orbital phase $ is given by v $ = v cos 9, where, 9 is the angle between the magnetic 
axis and line of sight, and is related to the inclination angle i, magnetic obliquity (3 and $ 
as 



cos9 = cosz cos (3 — smi sin/3 cos($ + 7r/2^ 



(2) 



Assuming % = 35° ± 5° and (3 = 58° ± 5° (IBrainerd fe Lamblll985l ). the observed modu- 
lation in the velocity of Fe XXVI line can be well explained by Equation [2] with the emitting 
region moving at a velocity of ~ 1100 km s^ 1 accreting on to a single pole. This velocity is 
close to the shock velocity of an accreting white dwarf of mass in the range of 0.6-0.7 M & . 
Thus, we interpret the observed variation in the amplitude of line shifts with phase as due 
to the aspect effect according to Equation [2J Similar phase depende nt line velocity varia- 
tion of few hundred km s _1 are reported in seve ral FUV lines too (IGansicke et al.l 11998 



Mauche fe Raymond! 1 19981 ; iHutchings et al.l |2002| ). Also, observations with Low Energy 
Transmission Grating onboard Chandra X-ray observatory showed that the emission line 
components are somewhat broader than the instrumental width, and the authors attribute 
this to the continu ously changing angle between the accretion column and the observer 
faurwitz et allbpOlh . 



The apparent lack of any velocity shift in the fluorescent line is in accordance with the 
standard emission models where this line is believed to be due to reflection of the hard X- 
rays from the surface of the white dwarf. The equivalent width of Fe fluores cent line (~ 150 
eV) c ombined with an N# smaller than 10 +21 cm -2 also supports this idea (jEzuka fe Ishida 
1999h . 
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In the above discussion we assume that the observed velocity shifts are real and not an 
artefact. This can be judged from Figure H] that clearly shows the varying line center of Fe 
XXVI and Fe XXV. Also the shifts in the line center values obtained from the overlapping 
phase bins (Figure E]) clearly show a sinusoidal trend modulated with the orbital period. 
The one artefact that might get introduced into our analysis is wrong identification of a 
particular line as the wavelength region around 1.6-2.0 A is filled with many DES and 
other lines. Though the resolution of Chandra HEG data does not allow us to fully resolve 
individual He-like lines of Fe, and wrong identification of a line is plausible, the effect of this 
is only to change the final velocity value but not the modulation. 



4. Conclusions 

We have presented here an analysis of high resolution Chandra observations of a pro- 
totype polar AM Her. The analysis of the average X-ray spectrum shows that it is well 
fitted by a multi-temperature, partially absorbed plasma emission APEC models. A mass 
of 1.15±0.05 M is derived based on the best fit of multi-temperature plasma model to the 
average X-ray spectrum of AM Her. The multi-temperature nature of the plasma is fur- 
ther confirmed by analysis of individual emission lines. The plasma diagnostics based on 
line-ratios of helium-triplets and Fe-L shell ions suggests a dense (> 2 x 10 12 cm -3 ) plasma. 
The helium-like triplets of O and Ne give a temperature of >2 MK, whereas the ratio of 
hydrogen-like line to helium-like line of Fe suggests a temperature of 12 keV. We have con- 
structed the DEM distribution of a AM Her, which is the first time that such an exercise 
has been carried out for a polar. The constructed DEM distribution using individual line 
fluxes shows that the DEM is continuous with a possible minimum near logT = 6.8. 

From the phase resolved spectroscopy, we have found a possible evidence for bulk motion 
of the ionized material in the accretion column of AM Her. The velocity of Fe XXVI line 
is modulated as a function of the orbital phase consistent with a single pole accretion. The 
Fe XXVI ions from the hottest plasma show a maximum velocity shift of ~1100 km s _1 , 
that is close to the shock velocity expected for a ~0.6 M white dwarf. The velocity shifts 
are observed to be much smaller for Fe XXV(r), and the hydrogen-like lines of S XVI and 
Mg XII. The variation in the velocity shifts of different ions is used to calculate the height of 
the accretion column from the white dwarf surface from where the line emission originates. 
An mass of ~1.0 M for the mass of the white dwarf is indicated by using the heights of 
Fe XXVI and Fe XXV line emitting regions estimated by the observed velocities of these 
lines. The multi-temperature mass estimate hence, gives an upper limit for the mass of the 
white dwarf owing to the bias introduced by using single temperature flow in the accretion 
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column. The density and temperature structure in the accretion column are also derived. 
No velocity shift is seen in the fluorescent Fe line, consistent with its origin due to reflection 
from the white dwarf surface. Higher resolution observations are required to detect the radial 
velocity modulation in this line. 

We would like to thank the referee for her/his comments helping in improving the paper. 
This research has made use of the Chandra Data Archive (CDA), part of the Chandra X- 
ray Observatory Center (CXC) is operated for NASA by the Smithsonian Astrophysical 
Observatory. VRR is pleased to acknowledge partial support from Kanwal Rekhi fellowship. 
The Gaussian convolution program used in the analysis is kindly supplied by David, P. 
Huenemoerder. 

REFERENCES 

Aizu, K., 1973, Prog. Theor, Phys., 49, 1184 

Arnaud, K., in ASP Conf. Ser. 101, Astronomical Data Software & Systems V, 1996, ed. 
G.H. Jacoby & J. Barnes (San Fransisco: ASP), pl7 

Brainerd, J. J., & Lamb, D. Q., 1985, in Cataclysmic Variables and Low-Mass X-ray Binaries, 
ed. D. Q. Lamb and J. Patterson (Dordrecht: D. Reidel), p247 

Bluementhal, G.R., Drake, G.W.F., Tucker, W.H., 1972, ApJ, 172, 205 

Burwitz, V., Reinsch, K., Barwig, H., 2001, Astronomische Gesellschaft Abstract Series, 18, 
192 

Canizares, C. R., Davis, J. E., Dewey, D., et al., 2005, PASP, 117, 1144 

Cropper M., Ramsay G., & Wu K., 1998, MNRAS, 293, 222 

Cropper, M., Wu, K., Ramsay, G., & Kocabiyik, A., 1999, MNRAS, 306, 684 

Dubau, J., Gabriel, A.H., Loulergue, M., Steenman-Clark, L., & Volonte, S., 1981, MNRAS, 
195, 705 

Ezuka, H., & Ishida, M., 1999, ApJs, 120, 277 
Gabriel, A.H., Jordan, C, 1969, MNRAS, 145, 241 

Gansicke, B. T., Beuermann, K., & D, de Martino, 1995, A&A, 303, 127 



- 15 - 

Gansicke, B. T., Hoard, D. W., Beuermann, K., Gansicke, B. T., Hoard, D. W., Beuermann, 
K., Sion, E. M., & Szkody, P., 1998, A&A, 338, 933 

Garmire, Gordon P., Bautz, Mark W., Ford, Peter G., Nousek, John A., & Ricker, George 
R. Jr., 2003, SPIE, 4851, 28 

Griffiths. N.W., & Jordan, C, 1998, ApJ, 497, 883 

Hellier, C, & Mukai, K., 2004, MNRAS, 352, 1037 

Hoggenwerf, R., Brickhouse, N. S., & Mauche, C. W., 2004, ApJ, 610, 411 

Hutchings, J. B., Fullerton, A. W., Cowley, A. P., k Schimidtke, P. C, 2002, AJ, 123, 2841 

Huenemoerder, David P., Canizares, Claude R., Drake, Jeremy J., & Sanz-Forcada, Jorge., 
2003, ApJ, 595, 1131 

Ishida, M., Matsuzaki, Keiichi., Fujimoto, R., Mukai, K., & Osborne, J. P., 1997, MNRAS, 
287, 651 

Kashyap, V. & Drake, J.J., 1998, ApJ, 503, 450 

Kashyap, V. & Drake, J.J., 2000, BASI, 28, 475 

Mauche, C. W., & Raymond, J. C, 1998, ApJ, 505, 869 

Mattei, J. A., Mayer, E. H., & Baldwin, M. E., 1980, S&T, 60, 180 

Mewe, R., Gronenschild, E.H.B.M., van den Oord, G.H.J. , 1985, A&AS, 62, 197 

Mouchet M., 1993, in White Dwarfs: Advances in Observation & Theory, ed. Barstow, M., 
p. 411 (Dordrecht: Kluwer) 

Mukai K., & Charles P. A., 1987, MNRAS, 226, 209 

Nauenberg, M., 1972, ApJ, 175, 417 

Oelgoetzty, J., & Pradhan, A.K., 2001, MNRAS, 327, L42 
Pike, CD., Philips, J.H., & Lang, J., 1996, ApJ, 464, 487 
Porquet, D., and Dubau, J., 2000, A&AS, 143, 495 

Saxton, C.J., Wu, K., Cropper, M., Ramsay, G., 2005, MNRAS, 360, 1091 



-16- 



Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Raymond, J. C, 2001 Spectroscopic 
Challenges of Photoionized Plasmas, ASP Conf. Series 247, pl59. ed. Gary Ferland 
and Daniel Wolf Savin 

Szkody, R, Raymond, J.C., Capps, R.W., 1982, ApJ, 257, 686 

Tapia, S., 1977, ApJ, 212, 125 

Terada, Y., Ishida, M., Makishima, K., Imanari, T., Fujimoto, R., Matsuzaki, K., & Kaneda, 
H., 2001, MNRAS, 328, 112 

Wesemael, F., Auer, L.H., Van Horn, H.M., Savedoff, M.P., 1980, ApJ Suppl, 43, 159 

Wu K., Changmugam G., & Shaviv G., 1995, ApJ, 455, 260 

Wu, K., Cropper, M., Ramsay, G., 2001, MNRAS, 327, 208 

Young, P., Donald P. S., & Stephen A. S., 1981, ApJ, 245, 1043 



This preprint was prepared with the AAS IATgX macros v5.2. 

Table 1. Spectral parameters for the best fit 4-T APEC model to the average spectrum of 
AM Her. The VEM is calculated for a distance of 91 pc. 



APEC 


kT 

(keV) 


log(T) 
(K) 


x 2 


Dof 


VEM 

(xl0 53 cm~ 3 ) 


1 

2 
3 
4 


1 3+°- 03 

u - lo -0.08 

u -°°-0.08 
9 7+O.6 
z - ' -0.4 

15-61I? 


6.i7l8:5§ 

6.8818$ 

7 ™+ - 08 

' • du -0.08 

8 26 +0 ' 12 
°- zu -o.io 


1747 
1716 
1654 
1624 


2767 
2765 
2763 
2761 


1 90+2.55 
2 c i +a38 

QR+5.19 
°- ou -1.26 

25.4l^ 8 
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Table 2. Measured fluxes, equivalent widths and Ionic column densities of the emission 
lines identified in the HETG spectrum of AM Her. 



Ion 


Wavelength 


Energy 


Flux 


Equivalent Width 




(A) 


(keV) 


(10 -5 photon cm~ 2 s _1 ) 


(eV) 


Fe XXVI 


1 701 +0.001 
L - ' oi -0.002 




9.3li| 


173 


Fe XXV (r) 


i ocro+0.001 

±.ooo_ 0001 


6.69ll8;gg 


7.6^1 


104 


Fel K a 


1 QQ7+0.002 
±.»oi_ 0001 


6 4n2+ - 006 


y -°-i.o 


151 


S XVI 


A 7 oc;+0.002 

4. <oo_ 0001 


2 fiiq+0-001 
z - oiy -o.ooi 


2.6181 


17.9 


S XV 


r n oo+0.001 

o.uoo_ 0016 


2 461 + - 001 


o.8i8i 


5.2 


Si XIV 


6-18418Z 


2 nn5 +0 - 001 

z - uuo -0.001 


9 c+0.3 

z -°-0.3 


12.2 


Mg XII 


8.424lg;»J 


1 /1 70+0.001 
-■-•^'^-o.ooi 


1 0+O.2 
1 - z -0.2 


3.9 


Fe XXIV 


1U.DZO_ 0001 


1 167+ - 001 
i - iD ' -0.001 


1 9+0.2 
1 - z -0.2 


3.9 


Fe XXIV 


11 181+0.084 


1 inq+ - 008 

1 - luy -o.ooi 


1 q+0.3 
± -°-0.4 


2.8 


Fe XXIV 


11 4Q9+ - 033 
11.4OZ_ 144 


1 ns^+ - 003 


o.8i8i 


1.7 


Fe XXII 2 


11.78 


1.055 


< 0.3 




Fe XXII 


11.92 


1.040 


< 0.6 




Ne X 1 


12 13Q+ - 055 

1Z.10»_Q QQ2 


1 n99+ - 005 
1 - uzz -o.ooi 


0+0.5 

2 - 8 -0.5 


5.5 


Ne IX (r) 


1 S 4.^9+0-136 

±0.40Z_ Q 0Q1 


q92+ - 009 
u - yzz -0.001 


1 5+ - 7 


2.5 


Ne IX (i) 


1 S tvM+0.225 


Q16+ 0015 
u.yio_ 0002 


o.8l8:5 


1.3 


Ne IX (f) 2 


1 S o44+ - 225 


0.90518:888 


< 0.6 




Fe XVII 


it; ni2 +0 - 022 


o.826i8:88^ 


1 o +0 - 7 


1.5 


Fe XIII 


16 02i +0 - 851 


n 774 +0.041 
u - ' '^-0.041 


1 9+0.9 
± - z -0.6 


2.2 


Fe XVII 


17 056 +0 - 065 


n 797+0.004 
u -' z ' -0.004 


1 8+ L1 


5.6 


Fe XVII 


17.100 


0.725 


6 +0 - 9 




VIII 


18.966t°;ig 


n 654+ ' 001 

U.DO4_ 0Q1 


9-8+11 


19.8 


VII (r) 


21.603±8:JS 


0.57418:88^ 


14.7U-S 


18.2 


VII (i) 


91 O1O+0.013 

zi.oiz_ 001 


o.569i8:881 


7 5+ 4 ' 2 
'•°-3.0 


9.6 


VII (f) 2 


22.064 


0.562 


<3.6 





1 Blend of Ne X and Fe XXII lines 

2 Flux is only an upper limit obtained by freezing the line center. 
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Table 3. Measured line velocity shifts of emission lines and corresponding heights of 

emitting regions of AM Her 



Line Id 


v (i) 

v avg 

(km s _1 ) 


y(2) 

<I>, const 

(km s _1 ) 


V (2) 
Y mod 

(km s _1 ) 


vf) 

(km s- 1 ) 


Height 
(km) 


Fe XXVI 


430±130 


220± 26 


790±40 


565±195 


932±110 


Fe XXV(r) 


850±90 


660± 70 


145±110 


770±75 


529±213 


s xvi 


510±80 


545±140 


110±120 


500±160 


315±220 


Mg XII 


1901^ 


240±160 


180±140 


280±195 


105±130 



x The error bars are based on Cash statistic with 90% confidence for a single parameter. 
2 Sinusoid+constant fit to the velocity shifts obtained from non-overlapping phase bins. 
3 A constant fit to the velocity shifts obtained from non-overlapping phase bins. 



Table 4. Parameters of shock in accretion column of AM Her 



M wd 




Vsh 


kT sfe 


h 


nf 




(M ) 


(Re) 


(km s _1 ) 


(keV) 


(km) 


(10 15 cm' 3 ) 


(km) 


0.39 


0.016 


765 


11.3 


394 


3.4 


2556 


0.50 


0.014 


922 


16.4 


691 


2.3 


2806 


0.60 


0.013 


1066 


21.9 


1067 


1.7 


3017 


0.75 


0.011 


1290 


32.1 


1889 


1.2 


3318 


0.91 


0.009 


1555 


46.7 


3312 


0.8 


3643 
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Fig. 1. — Chandra HETG spectrum of AM Her for the summed +1 and -1 orders. Only 
HEG spectrum is plotted for clarity. Data plotted are smoothed by a Gaussian convolution 
of width 0.0025 A HEG. The data are shown as histograms and the dashed line represents 
line-by-line fit to the data. Gaussians are used to model individual emission lines from 
various ionized elements and the continuum is modelled by an unabsorbed bremsstrahlung 
model. 



- 21 - 




log T 

Fig. 3. — VEM as a function of temperature. The dashed lines, solid lines and dash-dot lines 
corresponds to VEM of helium-like, hydrogen-like and Fe-L shell lines respectively. The thick 
solid line is the DEM derived using MCMC(M) reconstruction technique on a set of lines of 
hydrogen- like, helium-like of different atoms and highly ionized Fe line fluxes assuming solar 
photospheric abundances. The shaded areas are the 95% confidence limits corresponding to 
each temperature bin. 
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Fig. 4. — Gaussian fits (dashed lines) to the observed (histogram) Fe XXVI, Fe XXV (r,i & 
f) and fluorescent Fe lines. The orbital phase is shown in the inset. The width of each phase 
bin, A$ = 0.2. 
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Fig. 5. — Plot of the velocity shifts as a function of orbital phase for the emission lines 
Fe XXVI and Fe XXV as described in the text. The filled circles correspond to the velocity 
shifts observed in the non-overlapping phase intervals with A$ = 0.2. 
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Fig. 6. — Schematic representation of positions (in km) of different line emitting regions in 
the accretion column of AM Her. The emitting regions are labelled. The width (in km) of 
individual emitting region corresponds to the errors shown in Table [3j 



